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Structure, compressibility factor and dynamics of highly size-asymmetric binary 

hard-disk liquids 
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By using event-driven molecular dynamics simulation, we investigate effects of varying the area 
fraction of the smaller component on structure, compressibility factor and dynamics of the highly 
size-asymmetric binary hard-disk liquids. We find that the static pair correlations of the large 
disks are only weakly perturbed by adding small disks. The higher-order static correlations of 
the large disks, by contrast, can be strongly affected. Accordingly, the static correlation length 
deduced from the bond-orientation correlation functions first decreases significantly and then tends 
to reach a plateau as the area fraction of the small disks increases. The compressibility factor of 
the system first decreases and then increases upon increasing the area fraction of the small disks 
and separating different contributions to it allows to rationalize this non-monotonic phenomenon. 
Furthermore, adding small disks can influence dynamics of the system in quantitative and qualitative 
ways. For the large disks, the structural relaxation time increases monotonically with increasing the 
area fraction of the small disks at low and moderate area fractions of the large disks. In particular, 
"reentrant" behavior appears at sufficiently high area fractions of the large disks, strongly resembling 
the reentrant glass transition in short-ranged attractive colloids and the inverted glass transition 
in binary hard spheres with large size disparity. By tuning the area fraction of the small disks, 
relaxation process for the small disks shows concave-to-convex crossover and logarithmic decay 
behavior, as found in other binary mixtures with large size disparity. Moreover, diffusion of both 
species is suppressed by adding small disks. Long-time diffusion for the small disks shows power- 
law-like behavior at sufficiently high area fractions of the small disks, which implies precursors of a 
glass transition for the large disks and a localization transition for the small disks. Therefore, our 
results demonstrate the generic dynamic features in highly size-asymmetric binary mixtures. 

PACS numbers: 61.20.Ja, 61.20.Lc, 64.70.P- 



I. INTRODUCTION 

The physical mechanism of the liquid-glass transition 
remains unclear despite its fundamental and technologi- 
cal importance [lHJ]. Binary mixtures of particles, as a 
widely used model glass-former, play an important role 
in understanding the microscopic processes driving the 
glass formation during the past years. For example, the 
Kob- Andersen binary mixture has been used extensively 
in simulations to test theoretical predictions and analyze 
new phenomena within the glass community To 
study the glass transition, binary mixtures are often in- 
troduced merely as a means of suppressing crystallization 
since glass cannot be formed in one-component systems 
with conventional isotropic pair potentials. However, re- 
cent studies show that mixing two constituents with dif- 
ferent size not only offers a way to avoid crystallization 
but also possesses its own importance and exhibits very 
rich dynamical behavior yjj-[l9[. 

For moderate disparate mixtures, in which structure 
and dynamics of the two species are qualitatively the 
same, generic mixing effects close to the glass transi- 
tion have been identified [TlTfl4l | . Specifically, two qual- 
itatively different scenarios for the structural relaxation 
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have been demonstrated: Increasing the mixing percent- 
age of small particles leads to a speed up of long-time dy- 
namics for relatively large size disparity (say a size ratio 
5 = a s /<Ji = 0.6, where a s and oi indicate the diameters 
of small and large particles, respectively), while small 
disparity (e.g., 5 = 0.83) leads to a slowing down [ill]. 
Of special interest is the dynamic arrest in highly size- 
asymmetric binary mixtures, where the species' long- 
time dynamics can be quantitatively and qualitatively 
different. Moreno and Colmenero have demonstrated 
by molecular dynamics simulations fl7l Il8j that a bi- 
nary soft-sphere mixture with S — 0.4 shows anoma- 
lous dynamic features, including the sublinear behav- 
ior for mean-squared displacements, concave-to-convex 
crossover for the intermediate scattering functions by 
varying temperature or wave vector and logarithmic de- 
cay for specific wave vectors of density-density correla- 
tors. These striking dynamic features differ from the 
standard pictures of structural dynamic arrest in glass- 
forming liquids and resemble the mode-coupling theory 
(MCT) predictions for fluids confined in matrices with 
interconnected voids |2(J. They also suggest two mecha- 
nisms of the dynamic arrest for different types of parti- 
cles, i.e., dynamic arrest for the large particles originates 
from the competition between soft-sphere repulsion and 
depiction effects induced by neighboring small particles 
while for the small particles, it comes from the competi- 
tion between bulk-like dynamics (induced by neighboring 
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small particles) and confinement (induced by the slow 
matrix of large particles). Influence of the confinement 
effect on the glassy dynamics has also been discussed 
in a study of binary soft-disk mixtures with large size 
disparities [2l|. Voigtmann and Horbach [l5[ investi- 
gated the diffusion of a binary soft-sphere mixture with 
S = 0.35. Emergence of the anomalous and power-law- 
like diffusion for the small particles has been interpreted 
as a precursor of a double-transition scenario, which com- 
bines a glass transition and a separate small-particle lo- 
calization transition. This double-transition scenario was 
first identified using MCT for binary hard-sphere mix- 
tures [22|, HH . Very recently, it has been demonstrated 
by MCT calculations [l(| that multiple glasses can oc- 
cur even in the simplest size-asymmetric binary hard- 
sphere mixtures. These different kinds of glasses can be 
distinguished by considering whether small particles re- 
main mobile and whether small particles contribute sig- 
nificantly to perturbing the big-particle structure, and 
can be separated by sharp transitions that give rise to 
higher-order transition phenomena involving logarithmic 
decay of the relaxation process and power-law-likc be- 
havior of diffusion. In particular, the glass transition 
curve displays an "inverted shape" for sufficiently large 
size disparity. Finally, it should be noted that highly 
asymmetric binary mixtures also attract a lot of interest 
with the emphasis on th eir p hase behavior j23 - l28| and 
rheological properties [29l431j . 

In spite of the above interesting findings, the struc- 
ture, the compressibility factor, and in particular, the 
dynamics in highly size-asymmetric binary liquids are 
not well understood. In this work, we study highly size- 
asymmetric binary hard-disk liquids (6 = 0.2 and 0.15) 
by tuning the area fraction of the small disks and thus 
we can control the strength of the confinement effects for 
the small disks and hence the strength of the depiction 
effects for the large disks. We find that the static pair 
correlations of the large disks are only weakly perturbed 
by adding small disks. By contrast, the higher-order 
static correlations of the large disks, captured by the 
bond-orientation correlation functions, can be strongly 
affected. Accordingly, the static correlation length of 
the large disks deduced from the bond-orientation cor- 
relation functions first decreases significantly and then 
tends to reach a plateau as the area fraction of the 
small disks increases. The compressibility factor of the 
highly size-asymmetric binary hard-disk liquids shows 
non-monotonic change with the area fraction of the small 
disks, which can be rationalized by separating different 
contributions to it. We further find that adding small 
disks can influence dynamics of the system in quantita- 
tive and qualitative ways. For the large disks, the struc- 
tural relaxation time exhibits monotonic increase with 
increasing the area fraction of the small disks at low and 
moderate area fractions of the large disks and crosses 
over into "reentrant" behavior at sufficiently high area 
fractions of the large disks. This "reentrant" behav- 
ior strongly resembles the reentrant glass transition in 



short-ranged attractive colloids [32|, [33( (where two dis- 
tinct kinds of glasses exist, dominated respectively by re- 
pulsion and attraction) and the inverted glass transition 
in binary hard spheres with large size disparity 16j. By 
tuning the area fraction of the small disks, relaxation pro- 
cess for the small disks shows concave-to-convex crossover 
and logarithmic decay behavior, as found in other binary 
mixtures with large size disparity [l7l.[l8j. Moreover, dif- 
fusion of both species is suppressed by adding small disks. 
Long-time diffusion for the small disks shows power-law- 
like behavior at sufficiently large <f> s values, which implies 
that a glass transition for the large disks and a localiza- 
tion transition for the small disks can occur in our sys- 
tem. 

The article is organized as follows. In Sec. II the tech- 
nical details of our work is described. In Sec. Ill we first 
present the results of the glass transition in the absence 
of the small disks and then discuss the influence of adding 
small disks on the structure of the large disks, the com- 
pressibility factor of the system and the dynamics of the 
highly size-asymmetric binary hard-disk liquids. Finally, 
we conclude the paper in Sec. V. 



II. MODEL AND METHODS 

We perform event-driven molecular dynamics simula- 
tions (EDMD) [34( of quasi-binary hard disks with large 
size disparity under periodic boundary conditions. The 
number of the large disks is fixed to be Ni = 1000. To 
avoid crystallization, the diameters for the large disks are 
chosen equidistantly in the range 0.8 — 1.2 with an inter- 
val of 0.01, retaining the average diameter < 07 >= 1, 
th en the size polydisp ersity of the large disks is A = 
^J{< of > - < 07 > 2 )/ < (7/ >= 11.98%, where < • • • > 
is the average of the corresponding variable among all 
the large particles. We will show later that crystalliza- 
tion is indeed avoided and glassy dynamics develops on 
increasing density for hard-disk systems with such poly- 
dispersity. To model a binary mixture with sufficiently 
large size disparity, the size ratio <5 = a s / < ai > is cho- 
sen to be 0.2 (or 0.15). The area fraction of the large 
disks <\>i = nJ^jlitfj/^L 2 , where L denotes the sim- 
ulation box dimension, ranges from 0.76 to 0.784. By 
adjusting the small disk number N s , the area fraction 
of the small disks 4> s = ttN s <j 2 /4L 2 is varied from to 
0.06, then the largest numbers of the small disks are 1974 
and 3509 for S = 0.2 and 0.15 respectively. All the par- 
ticles have the same mass m, and the temperature T 
is irrelevant for athermal systems and we set it to be 
T = 1.0. Length, time and pressure are reported in units 
of < 07 >, y/m < 01 > 2 /ksT and ksT/ < 07 > 2 , where 
ks is the Boltzmann's constant. To properly create the 
initial configuration for the binary hard disks with large 
size disparity, the desired state with specific 0; and <j> s 
is generated by compressing a low-density liquid using 
the Lubachesky-Stillinger algorithm [35I437I ]. The start- 
ing configuration is then equilibrated for 10 5 before pro- 
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duction run for each state point and 8 independent runs 
are performed to improve the statistics. 



III. RESULTS AND DISCUSSION 

In this section, we first present the results of the glass 
transition in the absence of small disks. Then we discuss 
the effects of adding much smaller disks on the struc- 
ture of the large disks, the compressibility factor and the 
dynamics of the highly asymmetric binary hard-disk mix- 
tures. It is expected that the static pair correlations of 
the large disks will not be perturbed significantly by the 
small disks due to the large size disparity. However, we 
will show in the following that adding very small disks 
will strongly affect the higher-order static correlations of 
the large disks, the compressibility factor of the system 
and the dynamics of both species. 



Glass transition in the absence of small disks 
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FIG. 2: (a) The self-intermediate scattering function F l s (q Pl t) 
at varying <f>i for (f>s = 0. (b) <f>i dependence of a-relaxation 
time T a at cj> s =0. The red solid line is the result of the 
MCT power-law fitting r a ~ {<j> c — 0z)~ 7 with 7 = 3.8 and 
cf> c — 0.791. The green dashed line is the result of the Vogel- 
Fulcher-Tamman fitting r a ~ e"* 1 ' 1 * "*' 1 with D = 0.28 and 
d>c = 0.805. 
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FIG. 1: (a) The pair correlation function gi(r) and (b) the 
static structure factor Si(q) at low-g region at varying cj>i for 
0s =0. 



We first focus on the static structure of the liquids at 
4> s = 0, as characterized by the pair correlation func- 
tion and the static structure factor, which are defined 
as g a (r) = ^-^^-^—^ ^^5^ - \r jk \) and S a (q) = 
jf- < p a (q)pa(—(l) >, respectively. Here, a denotes the 
particle type (I for large disks and s for small disks), and 

Pa{q) = X^=i eiq ri wrtn Vj the position of particle j. 
The representative results are displayed in Fig. 1. It is 
seen that both gi(r) and Si(q) vary weakly when (f>i is 
increased from 0.76 to 0.784. However, the dynamics can 
slow down by orders of magnitude in the same density 
range, which will be shown later. Additionally, although 
crystallization is obviously prevented, the splitting sec- 
ond peaks in gi (r) and Si (q) become more apparent as (pi 
increases (this phenomenon is more apparent in the re- 
sults of Si (q) for our model) , suggesting the development 
of the locally preferred order on increasing density (i.e., 
the system is weakly frustrated and the local order is in- 
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deed hexatic in two dimensions). These static structural 
features on approaching the glass transition have been 
discussed extensively in the past few years [38[ . Evidence 
has been provided that the development of the locally 
preferred order can be connected to slow dynamics and 
dynamic heterogeneity of glass-forming liquids. 

The dynamic slowing down on approaching the 
glass transition can be characterized by the self- 
intermediate scattering function F"(q p ,t) = < 
£^ ie iq p -[ ri (t)-r 3 -(o)] >; where • = ^-J an d the wave 
number q p corresponds to the first peak of the static 
structure factor for the large disks. Fig. 2(a) shows 
the results for several densities. We observe that the 
relaxation slows down and becomes more stretched as 
4>i increases. In particular, a two-step decay emerges at 
sufficiently large </>/ values. This two-step process is a 
typical phenomenon in glass-forming liquids, and it re- 
flects the increasingly caged motion of particles at high 
area fractions. We define the structural relaxation time 
r a as F"(q P7 t = r a ) = 1/e. In Fig. 2(b) we plot 77 
as a function of 4>i. Clearly, t\ drastically increases as 
4>i increases. For the density range investigated, the 
simulation data can be well fitted by the MCT power- 
law: 77 ~ (4> c — 4>i)~ 7 , where <f> c is the MCT glass tran- 
sition point, or the Vogel-Fulchcr-Tamman (VFT) law: 
where D is the fragility index and </>o 
is the ideal glass-transition point. MCT fitting yields 
7 = 3.8 and <f> c = 0.791 and VFT fitting yields D = 0.28 
and 0o = 0.805. 

Therefore, we have demonstrated that the model pre- 
sented here possesses the ability to undergo the liquid- 
glass transition and reproduces essential aspects of a typ- 
ical glass-forming liquid. In the following sections, we 
will discuss the influence of the small disks on the struc- 
ture of the large disks, the compressibility factor and the 
dynamics of the binary mixtures. 
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FIG. 3: (a) The pair correlation function gi(r) and (b) the 
bond-orientation correlation function ga(r) at varying (f> s for 
8 = 0.2 and 61 = 0.784. 



B. Influence of small disks on the structure of the 
large disks 

Since the size disparity in this work is very large, it 
is expected that addition of the small disks will not per- 
turb the static correlations significantly. This is indeed 
seen in gi(r) (Fig. 3(a)) and Si(q) (data not shown) at 
varying <p s . Small but noticeable differences are also ex- 
hibited in Fig. 3(a): The first minimum of gi(r) slightly 
increases as <f> s increases and gi(r) on the whole changes 
mildly with <j> s . This observation is similar to the results 
of Ref. [IH for soft-sphere mixtures. Then it is interest- 
ing to assess whether the higher-order static correlations 
are strongly influenced by adding small disks since the 
importance of higher-order static correlations, in under- 
standing the glass transition, has been emphasized in re- 
cent years [39h42| . Because our system can form some 
locally preferred ordered structures (already indicated in 
Fig .1), we can use the bond-orientation correlation func- 
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FIG. 4: <f> s dependence of static correlation length ^ at i5 = 
0.2 and 4>i — 0.784. ^6 is obtained by fitting ge,{r)/g(r) with 
the OZ function (see text). The solid line is guide for the 
eyes. 
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acterize static correlations beyond the pair level. Here, 
ip J e is a sixfold bond-orientation order parameter, and de- 



fined as ipQ 



Ml 1 ,, 



where tia is the number 



of the nearest neighbors for particle j and determined 
by the Voronoi construction (34j, and 6^ is the angle 
between (r m — ry) and the x axis (particle m is a neigh- 
bor of particle j). These higher-order correlations re- 
veal more clearly the structural change of the large disks 
as <j) s increases, as shown in Fig. 3(b). The envelops 
of g\{r) I gi(r) can be well fitted by the Ornstein-Zernike 
(OZ) function r~ 1 ' 2 e~ r '^ 6 , and the obtained fitting pa- 
rameters ^6 stand for the static correlation lengths, which 
are shown in Fig. 4. It is found that first dramati- 
cally decreases as <j> s increases and then tends to reach a 
plateau at <fi s > 0.04. Therefore, although the addition 
of the small disks does not perturb the structure of the 
glass-forming liquids at the static pair level, the higher- 
order static correlations can be strongly influenced. Since 
we will further show that the dynamics of the large disks 
can be strongly affected by adding small disks, our results 
also imply the important role of higher-order static cor- 
relations in understanding the dynamics of glass-forming 
liquids. 



Compressibility factor of highly asymmetric 
binary hard disks 
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FIG. 5: Main: cj) s dependence of the compressibility factor Z 
at varying <j>i for 8 = 0.2. Inset: effect of 5 on Z at cj>i = 0.784. 

We now focus on the thermodynamics of the highly 
asymmetric binary hard-disk mixtures. For hard-particle 
systems, the only relevant thermodynamic quantity is the 
pressure P. The compressibility factor Z (and hence P) 
can be measured in an EDMD by total momentum ex- 
change from all the inter-particle collisions during a cer- 
tain time interval At, i.e., Z = N P k ^ T = l + X^At ~5$F5t 
with TV the total particle number. It is found that as <p s 
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FIG. 6: Different contributions to the total Z for 8 — 0.2 
at several <j>i values, s-s, s-1 and 1-1 indicate contributions 
from small-small, small-large and large-large particle colli- 
sions. The zoom in the inset highlights the contribution from 
the small-small particle collision. The results are similar for 
8 = 0.15. 



increases, P (data not shown) monotonically increases 
due to the increase of the total particle number (since 
Ni is fixed). However, the evolution of Z with (f> s is 
not monotonic, as shown in Fig. 5. As <j) s increases, 
Z first decreases and then increases, showing a minimum 
at some 4>™, which shifts to larger values as the size dis- 
parity enhances (0™ sa 0.01 for 8 = 0.2 and 0™ sa 0.02 
for 8 = 0.15), as evidenced in the inset of Fig. 5. Ad- 
ditionally, Z decreases at fixed tfii and 4> s values as 8 
decreases. 

To better understand the non-monotonic change of 
Z with 4> s , we can calculate different contributions to 
Z. Since in binary hard-particle systems, the follow- 
ing relation holds: Z = 1 + Z s — S + Z s -i + Zi_i, where 
Z s _ Sl Z s -i and Zi-i indicate contributions from small- 
small, small-large and large-large particle collisions and 

Z a _p — YliAt A 2JVTAt ° • ^ nc resur ts are shown in Fig. 6. 
We observe that as 4> s increases, first decreases sig- 
nificantly and then tends to reach a plateau at <j> s w 0.04, 
implying that the large-large particle collision frequency 
becomes low as (f> s increases and reaches nearly constant 
values at large 4> s values. Interestingly, the value of (f> s 
where Zi_\ approximately ceases to decrease is consis- 
tent with that where ^6 roughly ceases to decrease. By 
contrast, both Z s _i and Z s ^. s increase rapidly within the 
whole (j) s range. However, their absolute values are small 
as compared to due to the very large size disparity. 
Therefore, when <p s is small, the amount of the increase 
from Z s ^i and Z s ^ s cannot compensate the amount of 
the decrease from resulting into the first decrease of 
Z. As <j> s increases further, the decrease of Z\_i becomes 
mild. In this case, however, binary collisions occur more 
frequently between small and large particles and between 
small and small particles and the other two contributions 
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(Z s -i and Z s _ s ) still increase rapidly, which leads to the 
increase of Z at large </> s values. Moreover, as can be seen 
in the inset of Fig .6, the increase rate of Z s _ s becomes 
large as (f> s increases, implying enhancement of bulk-like 
effects for the small disks. We will show in the next that 
these effects can influence the dynamics of both species 
in quantitative and qualitative ways. 



D. Dynamics of highly asymmetric binary hard 
disks 
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FIG. 7: Particle trajectories during a time interval of t; at 
8 = 0.2 and <f> t = 0.784 for <j(> s = 0.01 ((a), (d)), S = 0.03 
((b), (e)) and cj> s = 0.05 ((c), (f)). (a), (b), (c) are for large 
disks and (d), (e), (f)) for small disks. 

We now turn to the results concerning the dynamics of 
the binary hard-disk liquids as <f> s gradually increases. To 
this end, we have studied the self-intermediate scattering 
functions and the mean squared displacements (MSD) 

< >= it < El"i M*) " r,(0)| 2 >, focusing on 

their^ s evolution. We note that Moreno and Colmen- 
ero [13, HH have studied the temperature dependence of 
the mean squared displacements and the density-density 
correlators in soft-sphere mixtures with several specific 
compositions. They have found several striking dynamic 
features (as described in the introduction part) at varying 
temperature or wave vector due to different mechanisms 
of dynamic arrest for both species. Here, we focus on the 
change of the dynamics as <p s increases. As mentioned 
above, the increase of </> s leads to the enhancement of 
depletion effects for the large disks and bulk-like effects 
for the small disks. This can be more clearly revealed 
in Fig. 7, where particle trajectories for both species 
are presented during a time interval of 77 (defined as 
F l s {q P ,t = T[) = 1/e) for three (f>i values. Obviously, the 
small disks can only move within the voids left by the 
large disks (comparing different species' trajectories in 
Fig. 7), i.e., they are under confinement. It is also seen 
in Fig. 7 that the bulk-like effects for the small disks 
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FIG. 8: <f>s evolution of the self-intermediate scattering func- 
tions of the large disks F^(q p ,t) at <5 = 0.2 for (a) <f>i = 0.77 
and (b) <f>i — 0.784. Here, q p corresponds to the first peak of 
Si(q) for <f> s — 0. The results are similar for S = 0.15. 



enhance as <f> s increases due to the increasing neighbor- 
ing small particles. In addition, it seems that relaxation 
processes for both species are not uniform in space and 
that the region where the large disks move more slowly 
corresponds to the place where there are stronger con- 
finement effects for the small disks. Therefore, we exp ect 
that some dynamic features observed in Ref . flH Il7l . Il8j 
will also emerge in our study as </> s increases. 

We first consider effects of adding small disks on the re- 
laxation process of the large disks. In Fig. 8, the change 
of the self-intermediate scattering functions with <j) s is 
presented for the large disks. We observe two different 
relaxation scenarios for small (or moderate) and large <pi 
values. For small 4>i values (e.g., 4>i = 0.77 in Fig. 8(a)), 
the time scales of F l s (q p ,t) for both initial and final de- 
cay increase monotonically as </> s increases. This means 
that the increase of the total area fractions at small <fii 
leads to a gradually slowing down of the dynamics for 
the large disks. When entering the regime of larger <f>i 
values, by contrast, another different relaxation scenario 
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FIG. 9: cj>a evolution of the self-intermediate scattering func- 
tions of the small diks F!!(q p , i) at 8 = 0.2 for (a) <f>i = 0.77 and 
(b) cj>i — 0.784. Here, q p corresponds to the first peak of Si(q) 
for <f) 3 = 0. The green dashed lines highlight the logarithmic 
decay of Fg(q p ,t). The results are similar for 5 = 0.15. 



occurs, as shown in Fig. 8(b). Specifically, when cf> s 
increases from to 0.06, the initial part of the struc- 
tural relaxation for the large disks slows down and the 
plateau value at intermediate times increases, but the 
time scale for the final decay first decreases and then 
increases, exhibiting "reentrant" behavior. As one start 
tracing the glass-transition line at the highest density, the 
non-monotonic (j) s dependence of 77 at (f>i = 0.784 means 
that the glass transition point for the large disks first 
becomes large and then shifts to smaller values upon in- 
creasing 4> s , which is reminiscent of the "reentrant" glass 
transition observed in colloid-polymer mixtures (short- 
ranged attractive colloids) [32, |33[ and the "inverted" 
glass-transition curve proposed recently in a MCT study 
of highly size-asymmetric binary hard spheres [l6j . Thus, 
a crossover for the structural relaxation of the large disks 
exists when the area fraction of the large disks is in- 
creased from small to large values. Turning to the re- 
laxation process of the small disks themselves, several 



FIG. 10: Main: <j> a evolution of relaxation times at <5 = 0.2 for 
(a) large and (b) small disks. Inset: Effect of S on relaxation 
times. 



interesting features are also found as <p s varies. It is seen 
in Fig. 9 that the shape of F*(q p ,t) is concave at small 
(j) s values and becomes more convex as <p s increases and 
that logarithmic decay also occurs for specific <f> s values. 
These dynamic features are strongly reminiscent of the 
findings in Ref. [13, 13- However, differing from the large 
disks, no crossover behavior occurs in the relaxation pro- 
cess of the small disks. To better characterize the </> s 
evolution of the relaxation process at different 0/ values 
for both species, we plot relaxation times for both species 
(defined as F"(q p ,t = T a ) = 1/e) as a function of <p s in 
Fig. 10. Clearly, r s increases monotonically at all 4>i val- 
ues as 4> s increases. For the large disks, the monotonic 
change of 77 with <p s at small <pi values crosses over into 
the non-monotonic behavior at large (f>i values. Moreover, 
we note that enhancing the size disparity (as 6 decreases) 
leads to a drop for both 77 and r s at fixed (f>i and 4>s val- 
ues and does not alter the monotonic change of t s and 
the crossover of ti as (j) s increases, as shown in the insets 
of Fig. 10. 

Thus, we have demonstrated that adding small disks 
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FIG. 11: 4>s evolution of the mean squared displacements at 
4>i = 0.784 and 8 = 0.2 for (a) large and (b) small disks. 



FIG. 12: The mean squared displacements of large (dotted 
lines) and small (solid lines) disks at varying (j>i and 4> s = 0.06 
for (a) 8 = 0.2 and (b) 8 = 0.15. 



can affect the relaxation process of both species in quan- 
titative and qualitative ways. <j) s evolution of the relax- 
ation process for the large disks is different for small and 
large 4>i values and several striking dynamic features for 
the small disks have been identified as <j) s varies. The 
influence of the small disks on the mean squared dis- 
placements can also be strong. In the first, we find from 
Fig. 11 that diffusion for both species slows down as <fi s 
increases, which is consistent with experimental results 
of Ref. [28( and implies that diffusion of both species is 
suppressed with increasing <p s . Moreover, we observe in 
Fig. 12 that long-time diffusion of the small disks dis- 
plays power-law-like behavior (e.g., < r 2 s (t) >oc with 
H < 1 at cj)i = 0.784 for t = 10 3 - 10 4 ) at sufficiently 
large (p s values, in good agreement with Ref. fl5| . and 
that the exponent y, decreases with increasing cf>i. This 
anomalous, power-law-like diffusion of small disks implies 
that a glass transition for the large disks and a localiza- 
tion transition for the small disks can occur at sufficiently 
large <p s values in our system. Furthermore, the power- 
law-like behavior becomes less evident for a smaller 8 
value (comparing Figs. 12(a) and 12(b)), indicating that 



the double-transition scenario occurs at larger <p s values 
as the size disparity enhances. 



IV. CONCLUSIONS 

In summary, we have investigated the structure, the 
compressibility factor and the dynamics of highly size- 
asymmetric binary hard-disk liquids. We have demon- 
strated that the addition of the small disks will not 
perturb the structure of the glass-forming liquids at 
the static pair level, but the higher-order static cor- 
relations can be strongly influenced, suggesting impor- 
tance of higher-order correlations in understanding dy- 
namic arrest of highly size-asymmetric binary mixtures. 
The compressibility factor of the system changes non- 
monotonically upon increasing the area fraction of the 
small disks and separating different contributions can ra- 
tionalize this phenomenon. Adding a much smaller com- 
ponent can influence dynamics of the system in quantita- 
tive and qualitative ways. For the large disks, the struc- 
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tural relaxation time exhibits monotonic change with the 
area fraction of the small disks at low and moderate area 
fractions of the large disks. "Reentrant" behavior for the 
relaxation of the large disks is displayed at sufficiently 
high area fractions of the large disks, which strongly re- 
sembles the reentrant glass transition in short-ranged at- 
tractive colloids and the "inverted" glass transition in 
binary hard spheres with large size disparity. By tuning 
the area fraction of the small disks, relaxation process 
for the small disks shows concave-to-convex crossover and 
logarithmic decay behavior, as found in other binary mix- 
tures. Diffusion of both species is suppressed by adding 
small disks, and in particular, the long-time diffusion for 
the small disks shows power-law-like behavior at suffi- 
ciently large cj> s values, which implies precursors of a glass 
transition for the large disks and a localization transition 



for the small disks. Our results can help to better under- 
stand dynamic arrest in highly size-asymmetric binary 
mixtures. 
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